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CN ! Abstract 

A waveguide coincides with a strip having two narrows of diameter e. Electron 
motion is described by the Helmholtz equation with Dirichlet boundary condition. The 
part of waveguide between the narrows plays the role of resonator and there can occur 
electron resonant tunneling. This phenomenon consists in the fact that, for an electron 
with energy E, the probability T{E) to pass from one part of the waveguide to the 
other part through the resonator has a sharp peak at E = Eres, where Eres denotes a 
Q^' "resonant" energy. In the present paper, we compare the asymptotics of Eres = -E'res(e) 

fH ■ and T{E) = T{E, e) as e — t- with the corresponding numerical results obtained by 

approximate computing the waveguide scattering matrix. We show that there exists 
a band of e where the asymptotics and numerical results are in close agreement. The 
numerical calculations become inefficient as e decreases; however, at such a condition 
the asymptotics remains reliable. On the other hand, the asymptotics gives way to 
^ I the numerical method as e increases; in fact, for wide narrows the resonant tunneling 

CD ' vanishes by itself. 

^ . Though, in the present paper, we consider only a 2D waveguide, the applicability 

I of the methods goes far beyond the above simplest model. In particular, the same 

ly-^ ' approach will work for asymptotic and numerical analysis of resonant tunneling in 3D 

O ' quantum waveguides. 



1 Introduction 

As an electron propagates in a quantum waveguide of variable cross-section, the waveguide 
narrows play the role of effective potential barriers for the longitudinal motion. The part 
of the waveguide between two narrows becomes a "resonator" , and there can arise resonant 
tunneling. It consists of the fact that, for an electron with energy E, the probability T{E) to 
pass from one part of the waveguide to the other through the resonator has a sharp peak at 
E = Eres, where Eres denotes a resonant energy. There are prospects for building a new class 
of nanosize electronics elements (transistors, electron energy monochromators, key devices) 
based on the phenomenon of resonant tunneling. To analyze their operation, it is important 
to know Eres, the height of the resonant peak, the behavior of T{E) for E close to Eres, etc. 

In [1], electron propagation was considered in a 3D waveguide with two cylindrical 
outlets to infinity and two narrows of small diameter ei and 62- The electron motion was 
described by the Helmholtz equation with Dirichlet boundary condition, radiation condition, 
and a wave number k between the first and the second thresholds. For the aforementioned 
characteristics of resonant tunneling, there were obtained asymptotics as 61,62 — ?■ 0. The 
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asymptotic formulas provide mainly a qualitative picture. In the present paper, we show 
that, being supplemented by some computations, the asymptotics can tell a useful quanti- 
tative information as well. Though the paper continues the studies in [1], nevertheless it is 
practically self-contained; let us explain its goal in detail. 

The asymptotic formulas in |lj include several unknown constant coefficients, which 
can be found by solving some boundary value problems independent of ei and 62- Here, in a 
model situation, we calculate approximately such coefficients, which enables us to take the 
asymptotics as numerical values of resonant tunneling characteristics for sufficiently small 
El and £2- This leads to the question which ei, 62 could be considered as "sufficiently 
small" ; in other words, where does the asymptotics work in a proper way? Though there 
is no universal answer for such a question, some examples give a grasp of what should be 
expected in analogous cases. To this end we calculate (also approximately) the scattering 
matrix and then compare the results obtained by the asymptotic and computational methods 
independently of one another. Generally, it can be predicted that numerical calculations will 
become inefficient as the narrow diameters decrease and the resonant peak turns out to be 
"too sharp" ; however, at such a condition the asymptotics should become more reliable. 
On the other hand, the asymptotics will give way to the numerical method as the narrow 
diameters increase; in fact, for wide narrows the resonant tunneling would vanish by itself. 
We observe these phenomena and show that there exists a band of the diameters, where the 
asymptotic and numerical approaches give compatible results. 

In the present paper, we consider a 2D waveguide that is a strip with two narrows of 
the same diameter e (see Fig. [2]). For 2D waveguides, the asymptotics of resonant tunneling 
characteristics (as £ — )■ 0) are published here for the first time; however, we do not prove the 
formulas in the paper. The reader could obtain the needed proofs by modifying arguments in 
[1] related to a 3D situation. Nonetheless, we analyze the structure of asymptotics in order 
to explain what constants in the asymptotic formulas have to be calculated numerically and 
how to do that by solving some boundary value problems independent of e. 

The paper consists of five sections. The mathematical model of the waveguide and 
statement of the problem are given in section 2. The asymptotic formulas are presented 
in section 3. Then, in section 4, we list the constants to be calculated in the asymptotics, 
describe the boundary value problems needed for the purpose, and present the methods for 
solving the problems numerically. In the same section we also describe a method we have 
used for approximate computation of the waveguide scattering matrix. Finally, section 5 
is devoted to comparing basic resonant tunneling characteristics obtained in two different 
ways, asymptotic and numerical, independent of one another. 

Though in the present paper we considered only a 2D waveguide, the applicability of 
the methods goes far beyond the above simplest model. In particular, the same approach 
will work for comparison asymptotic and numerical analysis of resonant tunneling in 3D 
quantum waveguides. 
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2 Statement of the problem 

To describe the domain G{e) in occupied by the waveguide, we first introduce two aux- 
ihary domains G and in M^. The domain G is the strip 

G = R X D = {{x,y) e R"^ : X e R = (-00, +00); y eD = {-1/2, 1/2)}. 

Let us define Q. Denote hj K a. double cone with vertex at the origin O that contains the 
X-axis and is symmetric about the coordinate axes. The set K (1 S^, where is a unit 
circle, consists of two simple arcs. Assume that Q contains the cone K and a neighborhood 
of its vertex; moreover, outside a large disk (centered at the origin) Q coincides with K. The 
boundary dfl of Q is supposed to be smooth (see Fig. [T]). 



Figure 1: The set f2. 

We now turn to the waveguide G{e). Denote by ^{s) the domain obtained from Q by 
the contraction with center at O and coefficient e. In other words, {x,y) G ^{s) if and only 
if {x/e,y/e) e fl. Let Kj and fij{e) stand for K and f2(e) shifted by the vector Vj = (x°,0), 
j = 1, 2. We assume that |x5 — X2I is sufficiently large so the distance from dKi fl dK2 to G 
is positive. We put 

Gie) = Gnni{e)nn2is) 

(see Fig. [2]). The wave function of a free electron of energy satisfies the boundary value 
problem 

Au{x,y) + eu{x,y) = 0, {x,y)eG{e), (2.1) 
u{x,y) = 0, {x,y)edG{6). 

Moreover, u is subject to radiation conditions at infinity. To formulate the conditions we 
need the problem 



Aviy) + \My) = 0, yeD, 
v{-l/2) = v{l/2) = 0. 



(2.2) 



Figure 2: The waveguide G{e). 



The eigenvalues of this problem, where q = 1,2,... are called the thresholds; they form 
the sequence = (vrg//)^, q = 1,2,.... We suppose that k"^ in (12. ip is not a threshold. 
Given a real k, there exist finitely many linearly independent bounded wave functions. In 
the linear space spanned by such functions, a basis is formed by the wave functions subject 
to the radiation conditions 

M 

(2.3) 

Kj=l 

M 



UM+m{x, y) 



Y,SM-,m,{k)e-^^^^^,{y) + 0{(i 



5x\ 



X — >■ — CXO, 



M 



+ SM+m,M+jik) e'^^'^'^jiy) + 0(e-^^), x +oo. 



Here M is the number of the thresholds satisfying A^ < fc^; m = 1, 2, . . . , M; Um = \/k'^ — A^; 

is an eigenfunction of the problem (12. 2p that corresponds to the eigenvalue A^ and is 
chosen so that 



^Jljlvra sin Amy, m even. 



(2.4) 



'2/lvjnCos\my-, m odd. 

The function Uj{x,y) = e^^^^'^ j{y), j = 1,...,M, in the strip G is a wave incoming 
from — oo and outgoing to +oo, while UM+j{x, y) = e~^'^^^'^j{y), j = 1, . . . , M, is a wave 
going from +oo to — oo. The scattering matrix 

S = \\Smj\\m,j=l,...,2M 

is unitary. The values 

M M 



ITT'J I ' 



Tm ^ ^ I S-m, M+j I 
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are called the reflection and transition coefficients, relatively, for the wave 11^ incoming to 
G{e) from — oo, m = 1, . . . , M. (Similar definitions can be given for the wave UM+m coming 
from +00.) 

In the present work we will discuss only the case (tt//)^ < k"^ < (27r//)^, i.e., /c^ is 
between the ffist and the second thresholds. Then the scattering matrix is of size 2x2. We 
consider only the scattering of the wave incoming from — oo and denote the reflection and 
transition coefficients as 



The goal is to find a "resonant" value kr = k^ie) of the parameter k corresponding to the 
maximum of the transition coefficient, and to describe the behavior of T{k,e) for /c in a 
neighborhood of kr{e) as e — i- 0. 

3 Outline of the asymptotics 

When deriving an asymptotics of a wave function (i.e. solution of problem (12. ip ) as e — )• 0, 
we use the compound asymptotics method (the general theory of the method was exposed, 
e.g., in [2], [3]). To this end we introduce "limit" boundary value problems independent of 
the parameter e. Put G{0) = G H Ki n K2 (Fig. [3]); thus, G{0) consists of the three parts 
Gi, G2, and G^, where Gi and G^ are infinite domains while G2 is a bounded resonator. 



R = R{k,e) 



sn{k,e)\\ 



T = T{k,e) 



su{k,e)\'^. 



(2.5) 



G 



1 




G 
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Figure 3: The set G'(O). 



The problems 



Av{x, y) + k"^ 



'v{x,y) 
v{x,y) 



0, 
0, 



{x,y) G Gj, 
{x,y)edG, 



(3.1) 
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where j = 1,2,3 and dGj is the boundary of Gj, are called the first kind limit problems. 
Solutions Vi and ^3 are subject to some radiation conditions at infinity and all three functions 
^1, V2, V3 satisfy some conditions at the corner points. All of the conditions will be formulated 
as required. 

Let us turn to the domains Qi and ^2 (see Fig. [1]). Problems of the form 

Aw{^„Vj) = F{^„r^,) inQ,, (3.2) 
w{^j,rij) = ondQj, 

are called the second kind limit problems. We seek solutions of the problems satisfying 

^(0' Vj) = (pj^"^""') as Pj ^ 00; 

here {Cj,Vj) ^-^s rectangular coordinates with origin at the vertex Oj of Kj, pj being the 
distance from (^j, rjj) to Oj and Uj the opening oi Kj, j = 1,2. 

In the waveguide G{e), we consider the scattering of the wave U{x,y) = e*'^^^\l/i(?/) 
incoming from —00 (see fl2.4p ). The asymptotics of the wave function is the main technical 
result. Although rather cumbersome, it will lead to much more explicit characteristics of the 
process. The wave function admits the representation 

u{x, y; e) = Xi,eix, y)vi{x, y; e) + 

+Q{ri)wi{e''^xi,6'^yi; e) + X2,eix, y)v2{x, y; e) + (3.3) 
+e(r2)w2(£:"^X2, £"^2; £) + X3,eix, y)v3{x, y; e) + R{x, y; e). 

Let us explain the notation and the structure of this formula. When composing the formula, 
we first describe the behavior of the wave function to the right of the narrows, where the 
wave function can be approximated by a solution ^3 of the problem (13. ip in G^. The solution 
is subject to the radiation condition 

W3(x,2/;£) ~ si2(£)e^'^^^^i(?/) as x ^ +00, (3.4) 

the element Si2{e) of scattering matrix being yet unknown. Problem (13. ip does not contain 
e, nevertheless ^3 depends on the parameter because of Si2{s). By X3,e we denote a cut-off 
function defined by 

X3,6{x,y) = (1 - &{r2/e)) lGs{x,y), 

where r2 = y^xfT"^ and {x2, y2) are the coordinates of a point (x, y) in the system obtained 
by shifting the origin to the point O2; Igs is the indicator of G3 (equal to 1 in G3 and to 
outside G3); 0(p) is a smooth non- negative function on the half-axis ^ p < +00 that 
equals 1 as ^ p ^ 5 and vanishes as p ^ 25 [5 being a fixed small positive number). Thus 
X3,e is defined on the whole waveguide G{e) as well as the function X3,e^3 i^i (13.30 . 

Being substituted to (12. ip . the function X3,e^3 gives a discrepancy in the right-hand 
side of the Helmholtz equation; the discrepancy is supported near the second narrow (to the 
right of it). We compensate the principal part of the discrepancy by means of the second 
kind limit problem in the domain ^2. Namely, the discrepancy is rewritten into coordinates 
{^2,^2) in ^^2 and is taken as a right-hand side for the Laplace equation. The solution W2 of 
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the corresponding problem (I3.2p has to be rewritten into coordinates {x2,y2) and multiphed 
by a cut-off function. As a result, there arises the term Q{r2)w2{e~^X2, s~^y2; e) in fl3.3l) . 

Now we substitute the sum of two obtained terms into (12. ip . The principal part of 
the corresponding discrepancy is supported in G2 near the second narrow. We compensate 
it by solving the problem (13.10 in G2 and obtain the term X2,e{x, y)v2{x, y; e) with 

X2,e{x,y) = (1 - e(£-Vi) - e(£-V2)) icA^^y). 

Then in a similar way there arise 

e{ri)wi{e~^xi,e~'^yi;e) and Xi,ei{x,y)vi{x,y;e). 

At the last step, we find the function vi that satisfies both the limit problem (13. ip in Gi and 
the radiation condition 

v^{x, y; e) ~ Suie)a{e)e"'^^<^I^{y) + s,2{e)^3{e)e-''''^^^{y) 

as a; — )■ —00. The coefficients a, /3 and the entries Su, S12 of the scattering matrix turn out 
to be uniquely determined by a relation between a and /3 that assures compensation of the 
principal part of the discrepancy arising in the problem in Gi, and by requirements 

si2ie)a{£) = 1, si2ie)l3{e) = su{e). 

The remainder R{x, y; e) is small in comparison with the principal part of (13. 3 p as £ — t- 0. 

We specify (13. 3p provided k"^ varies in an interval containing a unique simple eigenvalue 
/eq of the problem (13. ip in G2. 

1. Introduce a special solution V3 of the problem (13. ip in G3 satisfying 

V3(x, y) ~ {r2'^^'^ + ar2^'^)<l>(v92) as r2 

(here and below, {rj,ipj) are polar coordinates with center at Oj, j = 1,2; ^{(p) = 
cos{7np / u)) , and 

V3{x,y) ~ Ae'"'^'''^i{y) as x +00. 

These conditions define V3 uniquely. The constants a, A (depending on k and on the 
geometry of G3) have to be calculated. We have 

V3{x,y;e) = -^v-i{x,y). 

2. Consider a solution Wr of the homogeneous problem (13.20 satisfying 

,„ ^^-S (P"^" + (^P-'^^^M^) + 0(p-3-/-), as p -> 00, e > 0; 

- I ^^-V-$(7r-<^) + 0(p-3-/-), asp^oo,e<0. ^ ^ ^ 

The constants a, (3 (depending on VL) have to be calculated. One can prove that /3 7^ 
(cf. [1], proof of Proposition 3.4). We put 

w-(e,r/) = ^ («;.(e,r/) - Cr(e, r/)(p'^/" + «p-'^/")$(^) - r/)/3p-^/'^$(vr - ^)) , 
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where Cr is a cut-off function equal to 1 — 9 as ^ > and as ^ < 0, rf) = (r{—^, v)- 
We also put 



Then 



3. Remind that is a simple eigenvalue. Let Vq be an eigenfunction corresponding to fcg 
and normalized by j^^ {vqI"^ dxdy = 1. We have 



6ir^^'^$((/9i), as ri 0; 
62r2^'^$(7r - V92), as r2 0. 



^o(a;,y) ~ <( , _ , _ , ^ (3-6) 



In what follows we assume that 61 7^ 0; such an assumption is fulfilled, for example, 
if fco is the first eigenvalue of problem (13. ip . Since G2 is invariant with respect to the 
transformation {x, y) ^ {d — x, y), while d = Ix? — x^l is the distance between Oi and 
O2, one can prove that q := &2/&1 = ±1- Introduce special solutions V21, V22 of the 
problem f l3.ip in G2 satisfying 



((P - kDr-""^^ + cir[/") $((^1), as r, -> 0; 

C2r2^'^<l>(7r - (/}2), as r2 0, 

Vr^" + rfir^/") $(^1), as n ^ 0; 

-6ir2"''^'^ + (i2r2^'^ ) $(7r - V92), as r2 -)■ 0. 



V2i(a;,2/) 
^22{x,y) 



The coefficients Ci, C2, di, c/2 depend on and on G2- One can prove that Cjiko) = bibj. 
Then 

V2{x,y;e) = ^^^^{{bilie) - c?2^(^))v2i(x, + C2(5(e)v22(x, y)), 

C'lC2 



where 



7(e) = (e-^-/- - aa) , ^(5) = (« + - «')^'^/^) ■ (3-7) 

4. We have 

wi{^um;e) = {(Asn{e) + A)e-^/'^w-{^,,r],) + {aAsu{e) + aA)e^/^w+{^,,r],)) , 

where su is defined by fl3.8p below. 

5. Introduce a special solution vi of the problem (13. ip in Gi by vi{xi, yi) = w^{d — xi,yi). 
Then 

vi{x,y; e) = ^^^\i{x,y) + ^\i{x,y), 
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where 

s^,ie) = i2tb,C2r\ik^-kl)h\^ie)\^-iik^-kl)d2-b2C2)W)Sie) 



+hca{e)6ie) - (cirfa - C2d,)\5i6)\^)suie), (3.8) 
suie) = 2ihc2{-{e-kl)bMef + {{e-kl)d2-biCi-b2C2)i{e)5{e) 

+ {cid2-C2dr)5{efY\ (3.9) 

while 7(e), 5{e) are defined by (13. 7p . One can verify that jsnp + |si2p = 1 (cf. pLj). 

Analysis of (13. 9p shows that T = T{k,e) = |si2p has a sharp peak at /c = fcres, 

k^,,, = kl- 2a6?£2"/" + O(£'"/"+0, (3.10) 

where r = min{7r/u;,2 — a}, a being a small positive number. Suppose that k varies in a 
small neighborhood / of Kes-, I = {k : \k - Kes] ^ ce^'^/'^+P}, p > 0. Then (13.91) takes the 
form 

q{A{ko )/\A{ko) \f 
1 — ii^ 



.i2(^,e) = '^^"^T.'r^ ,2 \ (1 + o(^^)), 



"res 



r47r/a; 



where P = {2bil3'\A{ko)\')- . Hence 



n^,^) = ,2 (1 + 0(e^)). (3.11) 



1 + P2 



res 



e 



47r/aj 



The width of the peak at its half-height (the so-called a resonator quality factor) is 

T{e) = (3.12) 

4 Problems and methods for numerical analysis 

The principal parts of asymptotic formulas (13.101) - (13.121) for the main characteristics of 
resonant tunneling contain the constants 61, \A\, a, /3. To find the constants we have to 
solve numerically several boundary value problems. In this section, we state the problems 
and describe a way to solve them. We also outline a method for computing the waveguide 
scattering matrix S. 

To find bi, we solve the spectral problem (13. ip in G2 by FEM as usual. Let Vq be an 
eigenfunction corresponding to k^ and normalized by J^^ \ VQ{x,y)\'^ dxdy = 1 . Then bi in 
(13. 6p can be defined by 

$(0) ^ ^ Vo{e,0). 

Let us calculate \A\. In order to avoid dealing with vi, which increases at Oi, we 
introduce v = (vi — yi)/A, 

ar^^'^$((y9i) as r2 — )■ 0; 



+ ] ^^/^{y-^) + 0(e~*l^il) as xi -00, 
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where a = 2ilma/A. According to Lemma 4.1 in [1], Ima = so a = 2iA. Thus, it 
suffices to calculate a. Denote the truncated domain Gi fl {{xi,yi) : Xi > —R} by Gf and 
the artificial part of the boundary dG^ fl {{xi,yi) : xi = —R} by T^. Let V he a. solution 
of the problem 

AVixuyi) + k^V{x,,y,) =0, {x,,y,)eGf; 

V{x,,y,) =0, {xuyi)edGf\T^; (4.2) 

dnVixuyi) + mV{xi,yi) = 2iuie'''^^^i{yi), (xi, yi) G T^. 

We find V with FEM and put 

a= v^£-"/'^V^(-e,0). 

Pass to description of a boundary value problem for calculating a, /3 in (13. 5p . Denote 
Q n {{r,(p) : r < R} by and dfl n {{r,(f) : r = R} by T^. Consider the problem 

Aw{^,7i) =0, (e,^)efi''; 
w(e,r/) =0, (e,?7)e9fiV''; (4.3) 

If is a solution and C > 0, then 

||^;L2(r«)|Kr'||^?;L2(r^)||. (4.4) 
Indeed, substitute u = v = w to the Green formula 

= {dnU, v)QQR\r^ + {dnU + Cu, v)rR - Ciu, v)rR - (Vm, Vv)qr 

= {g, w)rn - C\\w; HT^W - \\Vw; L,{n'')f. 
From this and the obvious chain of inequalities 

C|k;L2(r«)f ^ Clk;^2(r«)f + \\Vw; L.in^W = {g,w)rR ^ lk;L2(r^)|| \\g;HT^)\\ 

we obtain (14.41) . Denote the left part of F'^ by F^ and the right part of F'^ by F;^. Let 
W be the solution of (iJD as C = n/uR, qI^^r = 0, gl^^R = {2t: /uj)R^^''^^~^^ip). Since 
the asymptotics (13.51) can be differentiated, Wr-W satisfies (gj]) with g = 0{R-^^^/'^'^-^). 
According to (14. 4p . 

\\Wr - W- L2(F^)|| ^ c— = c'i?-3-/- 

as i? +00. We find W with FEM and take 

/3 = "^^^^R""'^ = V^W{-R,0)R-/-. 
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Obviously, \\{wr - - (W - i?'^/'^<l>(^)); L2(r^)|| ^ c'R^'^^/^, therefore we put 

Finally, we outline the method of calculating the scattering matrix. Introduce the 
notation 

G{e, R) = G{€) r]{{x,y):-R<x<d + R}, 

rf = dG{e, R)n{{x,y) ■.x = -R}, = dG{e, R) D {{x,y) : x = d + R} 

for large R. We search the row (s^i, . . . , Sm,2M) of the scattering matrix s = s{k) defined 
by (12.31) . m = 1, . . . , M. As approximation to the row we take the minimizer of a quadratic 
functional. To construct such a functional we consider the problem 

AX^ + k^X^ = in G{e,R), 



= on i?) \ (rf u 



2 

M 

R 



M 

{dn + iQX^ = J2 ^{^, + C)e^^^("+^)v^,(y) on (4.5) 
i=i 

where C ^ \ {0} is an arbitrary fixed number, and ai, . . . , om are complex numbers. The 
solution Um to the homogeneous problem (12.11) satisfies the first two equations (14.51) . The 
asymptotics (12.31) can be differentiated so Um satisfies the last two equations in (14. 5 p up to 
an exponentially small discrepancy. As approximation for the row (s^i, • • • , Sm,2M) "we take 
the minimizer a^(R) — (a'j'(-R), . . . , a2A/(-^)) functional 

J>i,...,a2M) = ||A'^-e--%^-5^*^^a,e-^^vl/^.;L2(rf)f 

+ ll'^^-E,=i«*^+^-^^^^'^'''''^^^'^2(rf)ir, (4.6) 

where is a solution to problem (14.51) . As shown in [1], a°(i?, k) — Smj{k) with exponential 
rate as i? — )■ cxd and j = 1, . . . , 2M. To find the dependence of on ai, . . . , a2A/, we consider 
the problems 

Avf + k'^vf = in G{e,R), 

vf = on 9G(e,/2) \ (rf urf), 
idr, + tC)vf = ^(Tz/,+C)e^*"^%, onrf, 

{dn + tOvf = onrf; j = l,...,M; (4.7) 



and 



Au^^ + fc^^ = in G{e,R), 

vf = on 9G(e,i?) \ (rf urf), 

(9„ + ^CH^ = onrf, 

{dn + tC)vf = t{Tiyj + C)e^'"'^''^'^^'^j onrf; j = M + 1, . . . , 2M. (4.8) 
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Express by means of the solutions v"^ = v^j^ to problems fl4.7l) - fl4.8p . We have 
'^rriR + Ylij ^et US iutroduce the 2M x 2M-matrices with entries 



^3 ^r«' 



■^R 1 
2 



■^R 1 
2 



j, m = 1, . . . , M. We also put 

m=l,. . . ,M. The functional (14. 6 P can be written in the form 

J^{a, k) = {aS\k), a) + 2Re {T^{k), a) + g^{k), 

where J-"^ is the m-th row of the matrix J-"^ and (■, ■) is the inner product on C*^. The 
minimizer = a^{R, k) (a row) satisfies aP£^ + J-"^ = 0. Recall that we are searching the 
m-th row of the scattering matrix as m = 1, . . . , M. Along the same arguments one can 
prove that the found minimizer for m = M + 1, . . . , 2M serves as approximation to the m-th 
row of the scattering matrix. Therefore, as approximation s^{k) for the scattering matrix 
s{k) we take a solution to the equation s^S^ + = 0. 

When M = 1, i.e. fc^ is between the first and the second thresholds, we take ( = —vi- 
Then f f = f ^ = 0, £^ = (l/z/i)Id, and = —uiJ^^. 



5 Comparison of asymptotic and numerical results 

Let us compare the asymptotics fcres,a(^) resonant energy fc^es(^) approximate 
value fc^gg„(e) obtained by numerical method. Fig. H] shows good agreement with the values 
for 0.1 < £ < 0.5. We have 

l^res,a(^) ~ ^res,n(^) I /^res,a(^) < 10 ^ 

for 0.1 < e < 0.3 and only for e = 0.5 the ratio approaches 2.10"^. For e < 0.1 the numerical 
method is ill-conditioned. 

The difference between the asymptotic and numerical values is more significant for 
larger e because the asymptotics becomes not reliable. However, as the numerical method 
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Figure 4: Asymptotic description k'^es,ai^) (solid curve) and numerical description k^es,ni^) 
(dashed curve) for resonant energy k^^g{e). 





Figure 5: Transition coefficient for e = 0.2: asymptotic description Ta{k'^ — ^resa) (solid 
curve) and numerical description T„(A;^ — /c^es,n) (dashed curve) for transition coefficient 
r(/c^ — /c^es)- The width of resonant peak at height h: asymptotic A.a{h, e) = AA; numerical 
An{h,e)^BB. 
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Figure 6: The width A{h,e) of resonant peak for various h (dashed hne for numerical 
description, sohd hne for asymptotic description): hne 1 for h = 0.2; hne 2 for h = 0.5; hne 
3 for /i = 0.7. 
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Figure 7: Ratio An{h, e) / Aa{h, e) as function in e. The ratio is independent of h within the 
accuracy of the analysis. 
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shows, for e > 0.5 the resonant peak turns out to be so wide that the resonant tunnehng 
phenomenon dies out by itself. 

The forms of "asymptotic" and "numerical" resonant peaks are almost the same (see 
Fig. [5]). The difference between the peaks is quantitatively depicted in Fig. [61 Moreover, it 
turns out that the ratio of the width An{h,e) of numerical peak at height h to Aa{h,e) of 
asymptotic peak is independent of h. The ratio as function in e is displayed in Fig. [71 

Note that for e = 0.1, i.e., at the left end of the band where the numerical and 
asymptotic results can be compared, the disparity of the results is more significant for the 
width of resonant peak than that for the resonant energy. 
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